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Abstract 



The extremely small probability of tunneling through an almost classical 
potential barrier may become not small under the action of the specially 
adapted non-stationary signal which selects the certain particle energy Er. 
For particle energies close to this value, the tunneling rate is not small during 
a finite interval of time and has a very sharp peak at the energy Eft. After 
entering inside the barrier, the particle emits electromagnetic quanta and 
exits the barrier with a lower energy. The signal amplitude can be much less 
compared to the field of the static barrier. This phenomenon can be called 
the Euclidean resonance since the under-barrier motion occurs in imaginary 
time. The resonance may stimulate chemical and biochemical reactions in a 
selective way by adapting the signal to a certain particular chemical bond. 
The resonance may be used in search of the soft alpha-decay for which a 
conventional observation is impossible due to an extremely small decay rate. 



PACS numbers: 03.65.Sq, 42.50.Hz 
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I. INTRODUCTION 



Quantum tunneling under a potential barrier is a part of a variety of physical processes. 
The computation of the probability for a classically forbidden region has a certain peculiar- 
ity from the mathematical standpoint: there necessarily arises here the concept of motion 
in imaginary time or along a complex trajectory |]-|3[]. The famous semiclassical approach 
of Wentzel, Kramers, and Brillouin (WKB) |l[ for tunneling probability can be easily re- 
formulated in terms of classical trajectories in the complex time. The method of complex 
trajectories is also applicable to a non-stationary case [§,0]. This method has been further 
developed in papers [pHHl, where singularities of the trajectories in the complex plane were 



accounted for an arbitrary potential barrier (see also JTTJ] ) . Recent achievements in the 
semiclassical theory are presented in Refs. fli"2]-|iT)|l . 

A control of tunneling processes by external signal is a part of the field called quantum 
control which is actively developed now, see for example Ref. |T7|] and references therein. 



Excitations of molecules, when one should excite only particular chemical bonds |]18| , p0| , a 
control of acoustic waves in solids [19], formation of programmable atomic wave packets 



2T|, a control of electron states in heterostructures p2j, and a control of photo current in 



semiconductors |23|], are typical examples of control by laser pulses. 

Let us focus on main aspects of tunneling under nonstationary conditions. When the 
electric field £ cos fit acts on a tunneling particle of the initial energy E, it can absorb the 
quantum Q (with the probability proportional to the small parameter S 2 ) and tunnel after 
that in the more transparent part of the barrier with the higher energy E + Q. The pay 
in the absorption probability may be compensated by the probability gain in tunneling. In 
this case the system tends to absorb further quanta to low the total probability of barrier 
penetration. If fl is not big, the process of tunneling, with a simultaneous multi-quanta 
absorption, can be described in a semiclassical way by the method of classical trajectories in 
the complex time |]6|-[10|. When a tunneling particle of the energy E is acted by a short-time 
signal, the particle energy after escape is E + 5E, where the energy gain 5E = Nu, should 
be extremized with respect to the number of absorbed quanta N and the energy uj of each 
quantum. This extreme corresponds to the method of classical trajectories |M . 

According to the perturbation theory, the both probabilities, of absorption and emission 
of a quantum, are small being proportional to S 2 . After absorption, 5E is positive which 
enhances the total probability due to the increase of the tunneling rate. After emission, SE 
is negative and the particle should tunnel with a smaller energy in a less transparent part of 
the barrier; in this case there is no gain in the probability due to tunneling as for absorption. 
At first sight, tunneling in a nonstationary field cannot be assisted by an emission of quanta 
of this field since one should pay in probability twice. This conclusion is correct as soon as 
the nonstationary field is small and the perturbation theory is applicable. For tunneling the 
perturbation theory with respect to a non-stationary field stops to be valid already for very 
small fields |6|-|T0|,PH • What happens to emission-assisted tunneling when the nonstationary 
field is not very small? 

For a not very small nonstationary field the perturbation approach should be substituted 
by the semiclassical one. In this approach a multi-quanta absorption is accounted by classical 
trajectories in the complex time and the transition probability W ~ exp(— A) is expressed 
through the classical action A calculated by means of those trajectories. This differs from 
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the the stationary case described by the WKB formula W ~ exp[— Ao(E)]. The particle 
with the initial energy E emits quanta, "dives" on the level E— \ 5E |, and proceeds the 
penetration with this energy during the certain time interval 19. The "diving" process does 
not contribute much to the action A (in contrast to the perturbation theory) but the further 
motion adds Aq(E— | SE |) and the term 265E < 0. The last term is analogous to the 
occupation one, 5E/T, when the temperature T is substituted by 1/29 ||25|| . As shown in 
the paper, under some conditions of the nonstationary signal the "diving" process is instant 
in the complex time and weakly contributes to the total action. The condition of the instant 
"diving" is not necessary for the strong reduction of A. 

The emission-assisted tunneling exhibits a surprising feature: under action of the certain 
signal the tunneling rate can be not exponentially small since the action can tend to zero 
(A(Er) = 0) at some energy Er which depends on the signal parameters. The negative 
term 28SE in the action plays an important role in the dramatic reduction of the action 
A under certain conditions. According to this, the tunneling rate W ~ exp(— A) reaches a 
sharp maximum at the energy Er. The method used implies formally a small probability 
exp(— A) <C 1 but it can be of the order of 10 _1 instead of an extremely small value. In 
an anharmonic potential well a classical or quantum particle of the energy E, acted by the 
signal S cos fit, can exhibit a resonance if to adapt the frequency Q to the energy E |fZ6],|i| . 
Analogously, to get a maximum tunnel penetration at some given initial energy E of the 
particle, the external signal should be adapted to this energy. In the first case a motion occurs 
in real time but in the second case one should use imaginary (complex) time to describe 
the process. Traditionally, an action in imaginary time is called the Euclidean action which 
comes from relativity where the geometry is Euclidean in imaginary time. Extending this 
analogy, one can call the above tunnel phenomenon the Euclidean resonance. 

A scenario of the Euclidean resonance is the following: the particle with the energy E, 
before the barrier, leaves it having the exit energy E + 5E which is less than E; the position 
of the exit wave packet in time is close to the peak of the signal and the duration of the exit 
wave packet is roughly the under-barrier time 9; during this time interval the tunneling rate 
is not exponentially small for energies (Er — E)/Er ~ Aq 1 (Er). This means that under 
action of a specially adapted signals a particle can penetrate through very non-transparent 
barriers (big A ). One should emphasize, this process corresponds to a motion deep under 
a barrier but not a transition over its top. As shown in the paper, the signal amplitude, 
required for the Euclidean resonance, can be very small comparing to an electric field of the 
static barrier. 

In Sees. 1 1- IV the problem of the Euclidean resonance for a decay of the metastable 
state in the ^-potential well is solved exactly in the semiclassical sense. This means that 
for the wave function ip(x,t) = exp[iS(x,t)/H + iai(x,t) + iha2(x,t) + ...] (we use K in this 
relation) the action S, the first, and the second non-semiclassical correction <j\ and 02 have 
been found exactly in the analytical forms. The developed method enables to write down 
the exact expression for any a n . In Sec. V the semiclassical condition (a\ a 2 S> ...) is 
obtained. In Sec. VI the method of complex trajectories is shown to be valid for the problem. 
In Sees. VIII-IX the method of complex trajectories is applied to the Euclidean resonance 
in a smooth potential. In Sec. XI some particular manifestations of the resonance are 
considered. It can be used for a selective stimulation of chemical reactions. For alpha-decay 
of nuclei the experimentally registered energies of alpha-particles range in the approximate 
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interval of 2-10 MeV. For lower energies the tunneling probability is very low which makes 
observations to be impossible. However, there is no a principal ban for the soft alpha-decay, 
for example, with the energies of alpha-particles of the order of 1 KeV. The energy, counted 
from the bottom of the nuclear potential well, remains of the order of MeV. Adapting an 
external signal, with the duration of one or ten femtoseconds, to a possible energy level of 
an alpha-particle, one can try to search the soft alpha-decay. In principle, one can stimulate 
the conventional alpha-decay but this requires a very short signals of the order of 10~ 19 s. 



II. TRIANGULAR BARRIER 



As in Ref. f24|, we consider decay of the metastable state in the potential 

V(x) = V-Sq\x\- 



2(V-E) 



in 



5{x) 



(1) 



under action of a nonstationary electric field £(t). In the limit Sq — > 0, the energy E 
corresponds to the bound state in the ^-function potential well. The symmetric wave function 
ip(x, t) = exp[iS(x, t)+ia(x, t)] depends on the classical action S and the quantum correction 
a. The action S has the form |2l 



£2 rto rto 

S{x,t) dt x ir 00 - (t - tt) - J dt 2 h(t 

/to 
dtxh{t 



+ {V- E)t - Vt 



(2) 



where 



TOO 



^2m(V - E) 
So 



h(t) 



m 

So 



(3) 



and the function t(x,t) is defined by the equation 

mx . , . (t Q — t) 2 

— = i{t - t )r 00 + K0 J + 

Co I 



to 



dt^tx - t)h(ti) 



(4) 



For a symmetric pulse £(—t) = S(t) the moment t = plays an important role and it 
is convenient to introduce t (x) = —it o (x,0). According to Eq. (£|), the function t (x) is 
determined by the relation 



x 



So_ 

2m 



r-To 

2t t 00 - Tq - 2r 00 / d£h(£) 
Jo 



(5) 



By means of Eqs. (2) and (|5|) one can write down the first derivative of the action 

.ds(x, o) 



dx 



So[(p(r ) - Too] 



(6) 
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and the second derivative 



.d 2 S(x,0) 
i — 7-^ — = m 



1 + h{r ) 



dx 2 "" r 00 - r - T h(r ) 
Here r is expressed through x by Eq. (|5|) and 

V9(r) = r+ /"(WO 
./o 

In the static case, S(t) = 0, the action has the WKB form [|TJ 

2r 00 



(7) 



(8) 



-(y — e) 



(9) 



The point x = x±, where idS/dx = 0, is remarkable since at this point the particle comes 
out of the barrier. The exit point X\ is given by Eq. (|5|) with r = T\ where T\ is determined 
by the equation 

<p(tq) = r 00 (10) 
There is another remarkable point xi = xfa), where T2 can be found from the equation 

r 00 - r - r h(r ) = (11) 
The first derivative @ is finite at this point but id 2 S/dx 2 ~ {x% — x)^ 1 ^ 2 . 



III. THE PROBABILITY OF DECAY 



One can show, using Eq. (2), that 



n .0/ n s ^0 \.d 2 S(x,0) 
iS(x 1 ,t)=iS(x 1 ,0) + 



5m z 



dx 2 



Eq. ([TJD determines the duration of the exit wave packet 



f 2 

Sm 2 



,^ 2 g(x,0) 



-1/4 



(12) 



(13) 



Only those exit points X\ are physical for which the second derivative of the action is negative. 
The classical trajectory is given by the Newton equation 



m 



d 2 x(t) 
dt 2 



£ + £{t) 



One can show from the Hamilton- Jacoby equation [26|, that 



Im 



dS(x(t),t) m T „ ( dx\ 
dt 



1 lm \dt) +[^o + ^W]Imx(t) 



(14) 



(15) 
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The qualitative plot of ImS is shown in Fig. 1. We see that the imaginary part of the action 
is a constant at a real classical trajectory. One can find the probability J dx \ ip | 2 to find 
the particle outside the barrier, after an action of the pulse, by putting ip = exp(iS + ia), 



where a is a non-semiclassical part in the exact expression of the wave function fl24] . At 
every moment of time the x-integration is determined by the saddle point in the vicinity of 
the classical trajectory in Fig. 1 and does not depend on time, according to conservation of 
flux. For this reason, one can calculate the probability at t — when 



ip = exp 



+ ia 1 (x 1 , 0) 



X'l 



(16) 



According to Eq. (23) of Ref. fl2"4" |, the first non-semiclassical correction has the form 



i<7i(m, 0) = - in Ji - 2- [l + hin)} - ~ 

TOO ^ 



(17) 



Ignoring the finite duration 5t ([13]) of the exit wave packet and using Eq. (^) for the second 
derivative, one can obtain for the total probability to find the particle outside the barrier 
the following expression 



W(E, t) = t(V - E) exp(-A)) 



exp[-A(E)\ 



+ 4 

e £o v /(r 1 /i(r 1 ) + r 1 -r 00 )(l + Mn)) 







2m 



(x det - Xi) 



where Q(t) is the step function and 

A(E) = 2 ImSfa, 0) 

The first term corresponds to the conventional decay process (£(t) = 0) 0, when 



A = -(V-E)t 00 , 



It comes from the continuity equation 

dW 2 cty 

— — = 1m ip — — ; 

at m ox 



1> 



^m(V - E) 
p(x) 



exp ( i / pdx 



(18) 
(19) 

(20) 
(21) 



where the momentum has its classical form p(x) = ^j2m(E — V + Sqx). In Eq. (18) Xdet is 
the point of detection of the exit particle. Outside the barrier a correction of the trajectory 
(14]) by the non-stationary part S(t) is small. The first term in Eq. ([17]) is a short-time 
expansion of the full expression 1 — exp[— t(V — E) exp(— A )] corresponding to the multi- 
instanton approach. 



IV. THE NONSTATIONARY PULSE 



Now let us specify a particular shape of the pulse of the electric field 

8{t) = -£exp (-ftV - 2ft W) 



(22) 
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We consider further the limit Q 4 9 4 ^> 1 and, omitting the small width (Q 2 9) 1 ^ 9, one can 
write 



<p(r) = r - \6Q(t -9); A = ^ - exp (Q 4 9 4 ) (23) 

For the pulse fl2"2]) the function (p 2 (r) is plotted in Fig. 2. One can see that at some values 
of Tog three solutions of the equation (p 2 (r) = Tq are possible. According to this, there are 
different types of x-dependence of iS(x, 0) which are shown in Fig. 3. The thick solid curves 
in Fig. 2 represents the physical regime when the second derivative (^) is negative. The 
regimes, corresponding to the thin solid curves in Fig. 2, are not realized. This refers only 
to the vicinity of the moment t — 0; when a moment of time is not close to zero, the pulse is 
small and there is only one physical branch <^ 2 (r) = r 2 . We do not study here how the two 
regimes merge with time and consider only a time interval bigger than St when the equation 
(|17|) for the probability holds. 



V. THE SEMICLASSICAL CONDITION 

The semiclassical condition means a big action S comparing to the non-semiclassical part 



a As follows from the expression for a f24], it becomes big at the point x 2 where the 



second derivative of the action (^) is infinity. The first derivative (|6|) remains finite at this 
point which makes an essential difference with the static case (S(t) = 0) when the relations 
idS/dx = and id 2 S/dx 2 = oo are fulfilled at the same point giving rise to the Stokes lines 
f27j . In the present case there is no Stokes lines which makes the procedure of going around 
x 2 in the complex plane to be easy. This procedure enables to avoid the point X2, where the 
semiclassical condition breaks down, and to go over from one branch in Fig. 3 to another 
remaining within the semiclassical approach. For the pulse fl2"2"|), the semiclassical condition 
does not break down at the point x\. 

For the pulse Q2"2|), Eqs. fliTf ) and (^I|) result in the estimates 



1 n VlnSW 

1 n 2 9 ' 2 2tt 2 9 K ' 

Going around the point x 2 , one can put r = r 2 + 5r, where 5r =\ 5r \ exp(z0) and the 
phase <p varies between and tt. We have 

Too -r - T h(r ) ~ (r 00 - 9) 1 - exp (-4OWlnO 2 2 5rj (25) 

One can choose | 8r |~ (Q 2 9)~ l . In the semiclassical expansion a = o\ + + ■■■ it should 
be i > ffi > ff2 > - The estimations at the semi circle around T2, easily followed from 



24 



arc 



and 



ia(r 2 + St, 0) ~ \ In - J- (1 - A) (26) 

Z r 00 - V Zr 00 



■< ,n roo W9Hntt 2 9 2 

m{T > + 6T) ~(V^We r 00 - 9 (27) 
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Now the semiclassical condition o\ ^> <j 2 reads 

-^n¥lnn¥«(l/-^ (28) 

The method used implies that exp(— A) -C 1 but this parameter can be of the order of 10 _1 
instead of an extremely small value. 



VI. THE METHOD OF COMPLEX TRAJECTORIES 



The action A in Eq. ( |TTD can be written in the form 







dr 



(29) 



where x(tq) is given by Eq. (^). The relation (p9|) can be transformed into the following one 

f 2 



A = ^- 

m 



Jo 



(30) 



where T\ is determined by Eq. (fL0|). 

The result ([ft]) can be obtained by the method of trajectories in the imaginary (corre- 
sponding to the Euclidean norm) time t = ir 



A = 2 dr 



m ( dx\ . . _ . . ^ 



where 

V{x) = V - S x 

The trajectory satisfies the classical equation in the imaginary (Euclidean) time 

d 2 x{r) 



m 



dr 2 



with the boundary conditions 



dx(r) 



dr 



2(V-E) 



m 



-£o-£(r) 



x{n) = o 



dx(r) 
dr 



It follows from Eqs. (|33|) and (|3|) that 



dx 

m dr = ~ 



dm) 



The trajectory has the form 



x(t) 



2m 



t{ - t 2 + 2 / d£ drjh(rj) 
o Jo 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 
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The trajectory starts at the exit point x(0) = x%, with the energy E + 5E, and terminates 
at the enter point x(t{) = with the energy E. The action (fTO), in the method of classical 
trajectories, produces the same action ([!(]) followed from the Hamilton- Jacobi equation. 
The method of classical trajectories should be supplemented by the rule of selection of 
stable extrema 



d 2 S{x,0) 
dx 2 



< 



(37) 



j-i 



The energy of the exit particle E + SE can be found from the condition 



E + 5E = V - £ zi 



(38) 



and the energy gain is 



5E 



(V-E) 



TOO 



'00 



pti pro 

2 dr 
Jo Jo 



(39) 



VII. THE EUCLIDEAN RESONANCE 



The action fl3"0|) has the geometric interpretation 



A— — (Q + - Q_) (40) 
m 

where Q + and Q_ are the areas shown in Fig. 2 for the pulse (0). The plot of the action 
(fl0|) is shown in Fig. 4 where the thick curves correspond to the physical regimes under the 
condition fl37l) . The action satisfies the relation dA/dE = dA /dE, at t 00 = 6, where A is 



given by Eq. (20). The action A has to be inserted in Eq. (18) and at a given Too one should 
take the lower curve. In this sense the regions (b) and (c) in Fig. 4 are most important 
where 

A = 2(E R -E)9- E < E ext (41) 

Here the definitions are used 

e 2 s 2 A i 



E ext = E R --^-^l + — -Xj (A T - A) (A T - A) (42) 
where the resonance energy and the threshold pulse amplitude are introduced 

Within the exponential accuracy the probability (18) to find the particle outside the barrier 



is 



W ~ exp [-2{E R - E)6] 6 (E ext - E) (44) 
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As shown in Fig. 4, the expression fl44|) has the jump behavior as a function of E. At 
E ext < E the exit point x\ becomes complex resulting in the complex trajectory, defined by 
Eq. flnp in the real time, after exit the barrier. As follows from Eq. (0), the amplitude of 
the exit wave packet exponentially decays in time as for a particle with a complex energy. 
Hence the particle with E > E ext does not leave the barrier. 
The energy of the exit particle (for the lower curve in Fig. 4) 

E + 5E = V - 3(V - E R ) (45) 

turns to E at the upper end of the lower curve in Fig. 4. The duration of the exit wave 
packet fljjD in the case of the pulse fl22|) is 



St = 



S(V-E R ) 



1/4 



(46) 



Within the exponential accuracy the tunneling probability W ~ exp (— min {An ; A}) (18) is 
plotted as a function of energy in Fig. 5. At each amplitude of the pulse A there is a peak 



of W at the energy E ext (f42|). When A — > Ay the peak grows, reaches the value W ~ 1, 
and remains of the order of unity under further increase of the amplitude A > Ay as follows 
from Eq. (|44|). The function W(E) has a very sharp peak at E = Er which jumps down to 
W ~ exp(— A ) at E > E R . Formally, the action A becomes of the order of unity at E = E R 
and the semiclassical condition violates. This means that the energy E R — E should not be 
less than To get a classic resonance at some fixed energy in a non-harmonic potential 
one should adapt the frequency of the periodic signal to this energy and the oscillation 
amplitude depends on the strength of the signal. Analogously, to get a maximum tunnel 
penetration at some energy E R the external signal should be also adapted to this energy (by 
means of Eq. (fE|)) an d the intensity of the exit particles depends on the signal amplitude. 
In the first case a motion occurs in the real time but in the second case one should use the 
imaginary (complex) time to describe the process. Traditionally, an action in the imaginary 
time is called the Euclidean action which comes from the relativity where the geometry is 
Euclidean in the imaginary time. Extending this analogy, one can call the above tunnel 
phenomenon the Euclidean resonance. 

The non-stationary field fl22|) at Q A A 3> 1 can be called the instant signal since in the 
imaginary time it is well localized near the "moment" id. The Euclidean resonance can 
occur not only for instant signals, nevertheless, they have an advantage of enhancement of 
the effective amplitude S exp(f2 4 6* 4 ) ^> £ . For an instant signal, if the initial energy level E 
is fixed, the Euclidean resonance occurs at some resonance value Or defined by the condition 



e R = *?S^* ( 47) 



The effective amplitude of the signal A ( PB| ) should satisfy not a resonance but the threshold 
condition A > At- The probability of tunneling (f0|) can be rewritten in the form 

Q Q~ 

W~exp -3V3A (E)^- Q{0 R -O) (48) 
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where Aq(E) is the conventional action (p0|). 

The way of the particle under the barrier is shown in Fig. 6 in terms of its total energy. 
In the limit of big fi 4 # 4 the right-hand side of Eq. (^) has the form — £ + \£ 5(t — 6). 
According to this, the particle moves in the following way: (1) it starts at the point x = X\ 
(t = 0) and moves free, (2) it reaches the point x ~ approximately at the moment i9, (3) 
it gets the instant energy gain at the point x ~ at the moment i0, and (4) it arrives in the 
true point x = having the energy E. Since the non-stationary pulse is instant, the action 



(|3T|) does not change during the pulse and can be written as a sum of two parts 

A = A (E + 8E) + 298E (49) 

where the conventional action Aq is given by Eq. (|20|) under the condition Too = 0. Since 
5E < 0, the action increases Aq(E + 5E) > Aq(E), but the term —26 | 5E | provides a 
reduction of the total action which results in A = at E = Er. 



VIII. THE EUCLIDEAN RESONANCE FOR A SMOOTH POTENTIAL 



Suppose the static potential to have no singularities, like the potential ([!]). Such a 
potential is shown in Fig. 7(a). In this case the action A has the form 



A = -2 Im / dt 
'c 



m ( dx \ ' 
~2\di) 



V(x) + x£{t) + E 



(50) 



The contour of integration is shown in Fig. 7(b) where it starts at the point t = 0. The 
horizontal dashed line corresponds to the conventional case £(t) = 0. The classical trajectory 
satisfies the equation in the complex time 



m ^ + V\x) = 8{t) 
For a smooth potential let us choose the signal in the form 

te 



S(t) = -s 



t 2 + 6 2 



exp (-QH 2 ) 



(51) 



(52) 



We consider the case of the instant signal when Q 2 8 2 ^> 1 and S(t) is localized at the small 
part (2) of the order of (fl 2 8)~ 1 in Fig. 7(b). For the instant signal (^) the motion along 
the part (1) in Fig. 7(b) (the horizontal path in Fig. 7(a)) is free and x(t) = x(t), where the 
free trajectory x(t) satisfies the equation 



m ( dx\ 2 Tr ,_ N 

ai) +v{x) 



E + 5E 



(region 1) 



(53) 



A deviation of mdx/dt from mdx/dt starts at the region (2) in Fig. 7(b). If t belongs to 
the regions (2) (the vertical path in Fig. 7(a)) or (3) (the left horizontal path in Fig. 7(a)), 
one can write 



m 



dx 
dt 



* f dx 

dtxSit^+mi — 



(region 2) 



(54) 
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The classical momentum mdx/dt gets the jump J c dt£(t) after passing the part (2), but the 
coordinate x does not have a jump on the time scale (Q 2 9)~ 1 in the vicinity of the part (2). 
At the region (3) the motion again is free, x(t) = x(t), where the trajectory x(t) obeys the 
equation 

y(w) 2 + T/(£) = jE ( re S ion3 ) (55) 

and has the form x(t — ztq) (we omit a real time shift). Since, in the limit of an instant 
signal Q 2 9 2 3> 1, the function x(t) is almost continuous it follows that tq ~ 9. Otherwise, 
the real x(t) would jump going from (1) to (3). 
The classical energy gain 

6E = -fdt£(t)^ (56) 
can be calculated by inserting Eq. ([54D into Eq. (^). The result is 



I dt£(t)- — 
ieJc 2m 



1 2 



dt£{t) 



c 



(57) 



Since, close to the point id, exp(—Q 2 t 2 ) ~ exp[Q 2 6 2 — 2iQ 2 8(t—i6)} one can circle the contour 
C around the point t = i9, which produces 



n 2 e 2 \ I ® x \ 1 („ac^ 2 e 2 ^ 2 



The velocity dx/dr at r = 9 is real and negative. 

The duration (in the imaginary time) 9 of the part (1) in Fig. 7(b) can be expressed 
through 5E according to the classical formula 

9 = M f * - dX (59) 
V 2 J X0 yJV{x) -E-5E 

Eq. (|59D, together with the boundary conditions V(x{) = E+5E and V(xq) = E, determines 
the function SE(E, 9). At t — id the coordinate x ~ xq and Eq. (|!k3) results in 

It is remarkable that the value (60) coincides with the extreme of 5E in Eq. ( j58|) reached at 
£ = £ C (E, 9), where 

ir9£ ext (E, 9) exp (Vl 2 9 2 ) = y/-2m5E(E,9) (61) 

The action has the form 

A = 2V2 j 1 dxy/V(x)-E-5E(E,9) + 295E(E, 9); {£ = £ ext (E, 9)) 
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where the function SE(E,9) is determined by Eq. (|59|). The condition A = defines the 
resonance energy E R {6) and Eq. flBTl) gives the threshold amplitude <&r(#) = £ext(E R (9),9) 



£t{9) = P ^ } y/-2m5E(E R (9),e) (63) 

If the signal amplitude is fixed, the expression flp^) is valid only for one value of energy 
= i? ei . t which is determined by the equation £ = £ ext (E,9). If E ^ E ext the method of 
trajectories cannot be used and one should solve the Hamilton- Jacobi equation. £ ext is anal- 
ogous to the same parameter (fj2|) and corresponds to the peak in the escape probability W 
as in Fig. 5. One should note that whereas for a smooth potential the method of trajectories 
enables to obtain the result solely for one extreme energy, in the degenerated case of the 
singular potential ([I]) this method is applicable in an interval of energies. This is due to that 
for the singular potential (|l|) the maximum of W(E) is sharp and the analog of the dashed 
curve of Fig. 8 coincides with the solid curve in Fig. 5. This does not mean the method 
of trajectories to be useless for a smooth potential since the energy E ext corresponds to a 
maximum of the escape rate and Eq. (|62|) determines the dashed line in Fig. 8. In the case of 
a smooth potential it is impossible to obtain an analytical solution of the Hamilton- Jacobi 
equation to get the preexponential factor in the expression for the decay rate W as in the 
formula (18). 



IX. SMOOTH POTENTIAL WELL IN AN ELECTRIC FIELD 



Eqs. ( |5T?D and ( |5"2D allow to solve the problem of tunneling under action of an instant 
signal for a general semiclassical potential V(x). Let us consider a particular case of a 
smooth potential well in a weak electric field Eq. This means that in Fig. 7(a) xq <C x-y. In 
this situation one can use the approximation of triangular potential in Eq. ( (59| ) and ( 62|) 



49 

A = — [V - E - SE(E, 9)} + 295E(E, 



9= ^ m (V-E-8E) (64) 

On 



The extreme electric field ( |BTD is given by the formula 



7T 



S ext (E,9)exp (Q 2 9 2 ) 



2m(V -E) 



and the action ( p2|) has the form ([HP with the definition of E R (J4 



A = 29 

The extreme value E ext in Fig. 8 is 

Eext = E R 



{V - E) - 



MS 

6m 



(£ = S ext (E,9)) 



9 2 £ 2 



2m 



2 ( ix£ 2 fl 2 
3"U e 



G (St - £) 



(65) 



(66) 



(67) 
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where Er is given by Eq. (fl~3~D and the threshold value of the signal amplitude is 




If the energy level E in the well is fixed, the decay rate W is determined by Eq. ( pE8| ) with 
the resonance value of the signal parameter Or given by Eq. (|47|). 

The condition of applicability of the semiclassical approximation cannot be derived for a 
smooth potential with the all details, as the condition (p8|), since for this purpose one should 
know the solution of the Hamilton- Jacobi equation. Nevertheless, the typical feature of that 
condition (V — E)t ^> 1, where t is a characteristic time of the problem, can be used as an 
approximation. Since the typical time is (Q 2 9)~ 1 , the semiclassical condition reads 

n 2 e 2 < (v - E)e (69) 

The condition (^) does not contradict to a big value of Q 2 6 2 since in the semiclassical case 
(V — E)6 is a big parameter. 



X. OVER-BARRIER ESCAPE 



A metastable state in a well can decay also by the over-barrier transition picking up a 
number of quanta necessary to reach the top of the barrier (in our case there is no a thermal 
activation). This process is competitive with tunneling under the barrier. Let us estimate 
the probability of the over-barrier escape. The Fourier component of the signal fl52|) 



£u = —9£ exp 



oj 



tdt 



Att 2 / 7_ oc t 2 + e 2 



exp 



-Q 2 t 



ioj 
2Q 2 



(70) 



is determined by the saddle contribution (the horizontal path Imt = u/2Q 2 ) and by the pole 
integration at t = i9, if 9 < u/2Q 2 



4, = -11x68 exp (Q 2 6 2 - u6) - o) + 2iypH 



EujOVL 



oj 2 - (2{l 2 6) 2 



exp 



An 2 



(71) 



The probability of the over-barrier escape W Q -b is proportional to the absorption of N quanta 
which can be estimated as (£uj) 2N and can be written schematically in the form 

W . b ~ £ 2N jexp [2N(n 2 6 2 - lu6)] + exp (^-2N^j J (72) 

where iV = (V — E)/u and we ignore an interference of the two contributions. Since the 
typical signal amplitude is £ ~ exp(— Q 2 9 2 ), the probability (^|) reads 



W . 



exp \-26{V - E)\ + exp 



-2(V -E) 



n 2 e 2 



An 2 



(73) 
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The optimum frequency uo = 2Q 2 6 in the second term produces the same exponent as the 
first one. Therefore, the probability of the over-barrier escape, stimulated by the signal (p2|), 
can be estimated as 

W _ b ~ exp [-29(V - E)} (74) 

The exponent in the expression ([73]) is of the same magnitude as the action A far away from 
the resonance. For this reason, an escape mechanism under the action of the signal fl52|) is 
the tunneling process but not an over-barrier escape. From the mathematical point of view, 
it is understandable since the extreme of the action A is provided by the negative energy 
gain 5E but not by 5E — (V — E). 



XI. MANIFESTATIONS OF THE EUCLIDEAN RESONANCE 



In this section we point out some particular physical situations when specially adapted 
signals result or may result in a dramatically enhanced tunneling rate due to the Euclidean 
resonance. 



A. Selective ionization of atoms 



As a simplest example of the Euclidean resonance let us consider the ionization of the 



hydrogen atom by the signal (52). In this case one can use Fig. 7(a), considering the well to be 
extended in energy down to — oo. Suppose the constant electric field £ = 2 x 10 7 eV/cm. In 
this field the conventional tunneling rate is proportional to exp (— A ) = exp(— 76), where A 
is given by Eq. fl2"0| ) and (V— E) ~ 13.6eV. This provides the decay time (V—E) -1 exp (76) ~ 
10 9 years. Under the action of the signal with the amplitude £ > £t, where the threshold 
field is given by Eq. (|68|), the Euclidean resonance occurs at 



* fl = ^T^W 's (75) 



and the tunneling rate is determined by Eq. ( ^8|) which reads in this case 

W ~ lO^(e R -e)/e RQ ^ _ ^ (76) 



According to Eq. fl4"5D, (V — E — 5E) = 3(V — E). With the parameters chosen, Xi/x ~ 100 
and the duration of the exit wave packet ([46]) is 5t ~ 0.280 ~ 2.8 x 10 _15 s. Using the 
semiclassical condition (|6"§D, which now reads Q 2 6 2 200, one can take the signal parameter 
VL 2 2 = 15. This gives £t — 1.5eV/cm and Q^ 1 ~ 2.6 x 10~ 15 s. For a smaller constant 
electric field £ , the resonance value 6r and the threshold signal amplitude £t are also 
smaller. As follows from the relation ([TjJ), the uncertainty of the signal parameter 9 should 
be approximately 1% to be close to the resonance. When in the system there are a few types 
of atoms/molecules, one can make a selective ionization of a particular group of them. 
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B. Emission of electrons from a metallic surface 

Electrons can tunnel from a metallic surface when an electric field £q is applied perpen- 
dicular to it. In this case V is the work function and E can be at the Fermi energy or down. 
When the constant electric field is modulated by a successive set of signals (|52D separated 
by the time interval of the order of 9, there is the steady current of electrons distributed 
around the energy E + SE with the relative accuracy of approximately A$ . The current 
density can be estimated as j ~ envp/A Q , where n is the electron density in the metal and 
vp is the Fermi velocity. For the parameters of Sec. XI(A), the steady electric current from 
the metal surface is j ~ 10 8 A/cm 2 . If to use Eq. (^Dj) for A Q the current density can be 
written in the form j = aSo/e, where the conductivity is a ~ e 2 n/m(V — E). When the 
metallic sample is connected to a circuit, its electroneutrality recovers faster comparing to 
the time fl7|). 

C. Selective stimulation of chemical reactions 

A variety of chemical and biochemical reaction are characterized by the potential barrier 
for electrons which separates the electron state of the energy E before reaction and the state 
of the energy E\ after reaction (E x < E) |29|| . Usually, the tunneling rate for chemical reac- 
tions is negligible and they occur by thermal activation through the barrier which requires, 
in many cases, a reduction by a catalyst. Tunneling, stimulated by the Euclidean resonance, 
can provide new possibilities for reactions which are impossible under conventional condi- 
tions. Under the conditions of the resonance, the exit energy E + SE does not necessary 
coincides with E\ but during the motion of molecules towards each other the energy pa- 
rameters vary and at some distance the condition E + 5E = Ei can be encountered. Using 
the Euclidean resonance, one can select some particular bonds by the resonance to choose a 
type of reaction to be stimulated. 

By this method one can control also some nonvalent molecular conversions such as inter- 
nal rotations of a part of the molecule around a single bond, etc. In this case the Euclidean 
resonance stimulates tunneling of ions. For ionic tunneling the resonance Or is approximately 
two orders of magnitude bigger and St is smaller comparing to electron tunneling. 

D. Stimulation of alpha-decay 

In the problem of alpha-decay the potential energy of the alpha-particle and the nucleus 
is shown in Fig. 9, where the potential to the right of the nucleus has the form (3/x and the 
typical nucleus size is xq ~ lCT 13 cm. The conventional decay rate is W ~ exp(— Aq), where 

A = vr/^ (77) 

and m is the mass of the alpha-particle JT|. A properly adapted resonance signal can result 
in huge enhancement of alpha-decay. In a situation of the Euclidean resonance Eq. (|59|) 
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reads (with xq ~ 0) 



The action (E2) has the form 



4 (E + 5E) 3 / 2 ^ 



2/3 

A = ( STr/Jv^mfl ) - 20E (79) 



At the Euclidean resonance (A = 0), 



„ /27m £ + <5£ 1 

'« = W SEP ; T = 3 (80) 

where the energy E + 5E corresponds to the horizontal path in Fig. 9. The decay rate, 
within the exponential accuracy, is 

W ~ exp f-^Ao^-^-) 6 (6 R - 0) (81) 

The threshold amplitude of the signal (|63"D is 

For the decay Nd^ 4 -> Ce^ + a, one should take /? = 2 x 58e 2 and E = 1.9MeV. This 
gives 9r ~ 0.8 x lCT 19 s, A ~ 177, and (with fi 2 # 2 = 15) £t — 10 9 eV/cm. These parameters 
unlikely are suitable for an experiment. Nevertheless, the situation may be different. The 
energy of an alpha-particle E— Eq is always in the Mev-scale and the experimentally observed 
E varies approximately between 2 and 10 Mev for different nuclei. This does not exclude 
a situation when for some nuclei the energy of the alpha-particle E can be of the order of 
IKeV. Normally, alpha-decay of such nuclei hardly can be detected in a conventional way 
due to an extremely small decay rate since A Q ~ 1.3 x 10 4 flTTj). But under the conditions of 
the Euclidean resonance for E = IKev, the resonance value is 9r ~ 10~ 14 s (|78|) and one can 
obtain the flux of alpha-particles with the energy l/3KeV and the approximate duration 
10~ 15 s per one signal (^6|). One can take Vt 2 2 = 18 which gives St — 20eV/cm (|81|). The 
relative accuracy of the signal parameter near the resonance, according to Eq. (|8l"r) , should 
be 0.01% but 0, of the order of ten femtoseconds, relates with experimental facilities. In this 
way, one can search nuclei exhibiting the soft alpha-decay. 



XII. SIGNAL CONDITIONS 

The signal shape ( |5"2"D is not unique, nevertheless, we focus on this shape. The plot at 
Q 2 2 = 15 is shown in Fig. 10 where the vertical and the horizontal axes correspond to 
8{t)/E and t/0. A deviation of the signal from the dependence 

S(t) = -<4 exp (-ft 2 t 2 ) (83) 
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is less than 10%, however, the signal fl83|) does not lead to a Euclidean resonance, resulting 
only in a perturbative correction to the conventional formula W ~ exp(— A ). From math- 
ematical point of view, it is due to that the signal fl83|) does not have a pole at a complex 
time. We do not study here in full how fluctuations of the signal influence the resonance. 
For the resonance condition it should be {Or — 6) /9r < Aq 1 (^) which determines the signal 
accuracy with respect to the parameter 0. 

The amplitude of the threshold signal St is proportional to exp (— Q 2 2 ). Due to the 
semiclassical condition (|69"D, the parameter VL 2 6 2 <C Aq 1 can be big for very non-transparent 
(without a non- stationary signal) barriers, leading to a very small threshold signal. In the 
limit of a non-transparent barrier (A = oo), the amplitude of the threshold signal is S = 0. 
What is a restriction for this non-physical limit? 

Every physical system, strictly speaking, is not Hamiltonian and the energy levels always 
have a finite width due to interaction with other subsystems. In the classical limit the 
equation of motion contains the friction coefficient 7 

d 2 x dx 

m + mi dt + = ^ 

A physical origin of 7 depends on a particular situation. According to the theory of dissipa- 
tive quantum mechanics of Caldeira and Leggett [2l| (see also Ref. f28[), a friction disturbs 
the under-barrier motion resulting in the energy- loosing path (la) shown in Fig. 7(a). The 
frictionless path (1) and the friction path (la) correspond to the same time 9 and, due to 
this, the energy gain | SE | is less for the path (la). This makes impossible the resonance 
condition (A = 0) when the friction becomes sufficiently big 7 > 7 C . The approximate 
estimate of the critical friction is 7 C ~ 6~ x . Therefore, the Euclidean resonance exists for a 
dissipative system if the friction coefficient is not big 

7 < 7c ~ 0- 1 (85) 

The threshold amplitude can be estimated in a general case as St ~ £o ex P (—Q 2 2 ) where 
So is some typical electric field in the system. The restriction condition is 

n 2 6 2 < (V - E)6 < LzJl ( 86 ) 

7 



The main feature of the signal (|52| ) is its Lorentzian shape which results in the poles at the 
complex time. The Gaussian modulation is not obligatory and results only to the strong 
reduction of the threshold amplitude St- We do not consider here technical possibilities to 
tailor proper electromagnetic signals which should be in the approximate range of 1-100 fs. 
Recent achievements in formation of short signals are presented in Ref. . 



XIII. DISCUSSIONS AND CONCLUSIONS 

The probability of tunneling through an almost classical static barrier is exponentially 
small, according to the general principles of quantum mechanics. Nevertheless, a particle can 
penetrate through a very non-transparent barrier under the action of the specially adapted 
signal of an external field. The signal parameters select the particle energy Er. For particle 
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energies close to this resonance value, the tunneling rate is not small during some finite 
interval of time and has a very sharp peak at the resonance energy E R . The width of the 
peak (E R — E)/E R ~ Aq 1 (E r ) is connected to the static action A (E R ) which can be of the 
order of hundred or even bigger for a very non-transparent barrier. After entering inside the 
barrier the particle emits electromagnetic quanta and exit the barrier with a lower energy. 
It is reasonable to call this phenomenon the Euclidean resonance since the motion occurs in 
imaginary (complex) time. 

The Euclidean resonance can have various applications to the problems where quantum 
tunneling is an essential process. For example, chemical and biochemical reactions occur by 
thermal activation over a potential barrier which, in many cases, should be lowered by a 
catalyst. The potential chemical conversions, which cannot be realized due to high barriers, 
may get run by the effect of the Euclidean resonance. Since the resonance peak is very 
sharp, one can use the effect for selective stimulation of chemical reactions adapting the 
signal to some particular chemical bond(s). Another possible application of the Euclidean 
resonance is search of the soft alpha-decay when emitting alpha-particles have the energy in 
the interval of 0.1-10 KeV. The duration of the required signal is in the range 1-10 fs and 
the amplitude is of the order of 1-10 eV/cm. The soft alpha-decay, if it exists, cannot be 
observed in the conventional way due to extremely small decay rate. 

The phenomenon of the Euclidean resonance may occur not only in tunneling processes 
but also in the over-barrier reflection of particles in quantum mechanics [p]] and in classical 
reflection of macroscopic waves (electromagnetic, acoustic, etc.) from a spatially-smooth 
medium ||27|| . All these effects constitute a class of the Stokes phenomena since in a static 
case they have the common mathematic feature: the classical turning point generates the 



Stokes lines |27j where a jump of the sub-dominant (exponentially small) solution occurs. 
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FIGURES 




FIG. 1. The imaginary part of the action is a constant at the real classical trajectory x(t) after 
exit from the barrier. Only a part of the {x; t}-plane in the vicinity of the trajectory is shown. 
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FIG. 2. The plot of the function ip 2 (r) consists of the two parabolas t 2 and (r — \9) 2 connected 
by the vertical line. The solution r = t% of the equation (p 2 (r) = Tq determines the exit point x\ 
where dS(x,0)/dx = 0. Only the thick curves are physical. The same graph gives the geometrical 
condition (the equality of the areas Q + = Q ) for the Euclidean resonance. 
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x 2 x 
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(c) 



(d) 



FIG. 3. The action for different domains of Tq in Fig. 2. x\ corresponds to the stable regime 
with id 2 S/dx 2 < denoted in Fig. 2 by the thick solid curves. 
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W(E) 








E V E 



FIG. 5. At each amplitude of the signal, there is a peak of W at the energy E ext . When A > \r, 
the peak reaches its maximum value at E = Er. When E > Er, W is determined by the WKB 
formula. 
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V(x) 




FIG. 6. The total energy of the particle under the barrier. The particle starts at the point x\ 
(r = 0), reaches the point x ~ at the "moment" r ~ 9, acquires instantly the energy gain, and 
reaches the point x = having the energy E. At the horizontal path the particle energy is E + 5E. 
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FIG. 7. (a) The path of the particle under the barrier in terms of its total energy. The signal 
acts at the "diving" (vertical) part. The dashed curve shows the particle energy in the damped 
case, (b) In terms of the complex time, the part (1) represents the horizontal path in Fig. 7(b), 
the part (2) - the vertical path, and the part (3) - the left horizontal path. 
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FIG. 8. The tunneling rate W ~ exp(— A) at each fixed amplitude of the signal £ < 6t is a 
curve with the maximum at E = E ext . E ext = Er at £ = £t- W reaches its absolute maximum at 
E = Er under the condition £ > £t- 
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FIG. 10. The signal (52) is drawn by the solid curve at Q 2 9 2 = 15. The vertical 
and the horizontal axes correspond to £(t)/£ and t/9. The dashed curve shows the signal 

£{t)/£ = -(t/9) [i + (t/ef]-\ 



31 



